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Abstract. The non-commutative Central Limit Theorem (CLT) introduced by Speicher in 1992 
states that given a sequence (&i)i>i of non-commutative random variables satisfying the commuta- 
tion relation 

bib) = s(j, i) b) b\, s(j,i) 6 {-1,1}, 

the *-moments of the normalized partial sum Sn = (bi + . . . + 6jv)/V r /V are given by a Wick-type 
formula refined to count the number of crossings in the underlying pair-partitions. When coupled 
with explicit matrix models, the theorem yields random matrix models for creation and annihilation 
C$ operators on the q-Fock space of Bozejko and Speicher. 

In this paper, we derive a non-commutative CLT when the commutation relation is generalized 

£ — b\b) = n e > te (j,i)bj bi, fv >e (j,i) 6 R. 

The statistics of the limiting random variable are a second-parameter refinement of those above, 
jointly indexing the number of crossings and nestings in the underlying pair-partitions. Coupled 

t^J with analogous matrix constructions, the theorem yields random matrix models for creation and 

annihilation operators on the recently introduced (q, t)-Fock space. 

& 

■i— > 

3 1. Introduction 

In non- commutative probability, probabilistic interpretations of operator algebraic frameworks 
j>. give rise to non-commutative analogues of classical results in probability theory. The general 

setting is that of a non-commutative probability space (A, tp), formed by a *-algebra A, containing 
the non- commutative random variables, and a positive linear functional cp : A — > C, playing the 
role of expectation. A particularly rich non-commutative probabilistic framework is Voiculescu's 
free probability |10j , which has been found to both parallel and complement the classical theory 
(see in-depth treatments in [21 [HI [11] ) . Whereas free probability can be seen as characterized by the 
absence of commutative structure, a parallel - albeit somewhat slower - development has targeted 
non-commutative settings built around certain types of commutation relations. 
J> In [9 1, Speicher showed a non-commutative version of the classical Central Limit Theorem (CLT) 

for mixtures of commuting and anti-commuting elements. Speicher's CLT concerns a sequence of 
elements 6i , t»2 - • • £ A whose terms pair- wise satisfy the deformed commutation relation bibj = 
s(j,i)bjbi with s(j,i) S {—1, 1}- It is not a priori clear that the partial sums 

bi + . . . + b N 



should converge in some reasonable sense, nor that the limit should turn out to be a natural 
refinement of the Wick formula for classical Gaussians, but that indeed turns out to be the case. The 
following theorem is the "almost sure" version of the Central Limit Theorem of Speicher, presented 
as an amalgamation of Theorem 1 of [S] and Lemma 1 of [D]. Throughout this paper, ,5^2 (2n) 
will denote the collection of pair-partitions of [2re], with each 'f € ^2(2^) uniquely written as 
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'f = {(wi, zi), . . . , {w n , z n )} for w\ < . . . < w n and Wi < %i (i = 1, . . . , n). For further prerequisite 
definitions, the reader is referred to Section [2j 

Condition 1. Given a *-algebra and a state ip : .etf — > C, consider a sequence {6j}j g N of elements 
of stf satisfying the following: 

(1) for all i £ N, = ip(b*) = 0; 

(2) for all for all i, j £ N with i < j and e, e' £ {1, *}, <f{b\bf ) = ipipff-); 

(3) for all n G N and all . . . ,j(n) £ N, e(l), . . . , e(n) £ {!>*}; the corresponding mixed 
moment is uniformly bounded, viz. \<f(Y!7=l bjU))\ — ° n f or some non-negative real a n ; 

(4) ip factors over the naturally ordered products in {6i}jgN, in the sense of Definition^ 
Assume that for all i / j and all e, e' £ {1, *}, b\ and b e - satisfy the commutation relation 

brf = s(j,i)bpl s(j,i) £ {-1, 1}. (2) 



Theorem 1 (Non-commutative CLT [9]). Consider a non- commutative probability space (A, (p) 
and a sequence of elements {6i}igN in A satisfying Condition^ Fixing q £ [—1,1], let the com- 
mutation signs {s(i, j)}i<i<j be drawn from the collection of independent, identically distributed 
random variables taking values in {—1, 1} with E(s(i, j)) = q. Then, for almost every sign sequence 
{s(i, j)}i<i<j, the following holds: for every n £ N and all e(l), . . . , e(2n) £ {1, *}, 

lim rtsfK..sf n ~V) = 0, (3) 

TV— »oo 



lim <p(Sff> . . . Sf n) ) = Yl 9 Cr °" (r) II <P(b e{wi) b e ^) , (4) 

with Sn £ A as given in ([T]) ; "V = {(wi, Z\), . . . , (w n , z n )}, and where cross(y) denotes the number 
of crossings in (cf. Definition^. 

The moment expressions in © and @ may seem familiar. Indeed, this stochastic setting with 
the average value of the commutation coefficient set to q turns out to be, from the point of view 
of limiting distributions, equivalent to the setting of bounded linear operators on the the q-Fock 
space JP q (Jt?) of Bozejko and Speicher [3]. Specifically, given a real, separable Hilbert space Jif 
and two elements f,g£ Jf?, the creation and annihilation operators on ^ q ( r ^f), a q (f)* and a q (g) 
respectively, satisfy the q-commutation relation: 

a q(f) a q(sT ~ q a q (g)*a q {f) = (/, g) jf 1- (5) 

The mixed moments with respect to the vacuum expectation state <p q of these operators are given 
by a Wick-type formula which, compared against ^ and Q, yields that for a unit vector e in 
JT, limjv^oo v(S e i 1] . ■ ■ S^ n) ) = (p q (a q (e)<V . . . a q {ef^) for all n £ N and e(l), . . . , e{2n) £ {1, *}. 
As described in [9] , Theorem [I] can be used to provide a general asymptotic model for operators 
realizing the relation ([5]), thus providing non-constructive means of settling the conjecture in [BJ. 

Finally, any sequence {fri}ie[n] satisfying Condition [l] has a ^representation on s^ n := A42(M)® n , 
where A^BO denotes the algebra of 2 x 2 real matrices. Matricial models for operators satisfying 
the canonical anti-commutation relation, i.e. the fermionic case corresponding to 9 = — 1 in (pjb, 
are well known and are provided by the so-called Jordan- Wigner transform (see e.g. [5] for its 
appearance in a closely-related context). By extending the transform to the setting where there 
are both commuting and anti-commuting elements and by applying Theorem [T], Biane pQ obtained 
a random matrix model for operators satisfying the (/-commutation relation By replacing 

2x2 matrices with 4x4 block-diagonal matrices, Kemp [7] similarly obtained models for the 
corresponding complex family (a(f ) + ia{g)) /y/2. To describe the extended Jordan- Wigner model, 



we make the identification stf n = M2^(M) and let the * operation be the conjugate transpose on 
At2™(K)- Furthermore, let (p n : A^QR) — > C be the positive map a \-t (aeo, eo) n , where (, ) n is the 
usual inner product on W 1 and eo = (1, 0, . . . , 0) an element of the standard basis. 

Lemma 1 (Extended Jordan- Wigner Transform [lj). Fix q G [—1,1] and consider a sequence of 
commutation coefficients {s(i,j)}i<j drawn from {—1,1}. Consider the 2x2 matrices {o~ x } x ^, 7 
given as 





' 1 


" 




" 


1 " 


0~x = 





X 


7 = 









and, for i = 1, . . . , n, let the element b n ^ G A^2(C)® n be given by 



V» = °~s(l,i) ® °~s(2,i) 



0-1 • 



(6) 



Then, for every n G N, the non- commutative probablity space (A n ,ip n ) and the elements b ni i,b n} 2, 
■ ■ ■ , b n ,n G An satisfy Condition^ 

1.1. Main Results. This article derives a general form of the Non-commutative Central Limit 
Theorem of [9]. The setting now concerns a sequence {&i}jgN of non-commutative random variables 
satisfying the commutation relation 

b\b< = ^ >e (j,i)bpt with e,e' G{1,*}, /v, e (i,j)eR, (7) 

for i 7^ j. The consistency of the above commutation relation is ensured by requiring that for all 
i < j and e, e' G {1, *} 

Ml.ifoi) = T^v > Mi,*(*>i) = (8) 



(9) 



1 

Me,e'(«>i)' 

where i > is a fixed parameter that will appear explicitly in the limits of interest. The reader is 
referred to the beginning of Section [3] and the Remark [3] of Section [4] for a discussion of the relations 
((8j) — ([9]) ■ The conditions underlying the extended non-commutative CLT are now the following. 

Condition 2. Given a *-algebra and a state (p : srf — > C, consider a sequence {&i}jgN of elements 
of stf satisfying the following: 

(1) for all i G N, ^(64) = ip(b*) = </>(&A) = cp(b*b*) = ip(b*bi) = 0; 

(2) for all i,j G N and e, e' G {1, *}, p(6fbf ) = tp(V$); 

(3) for all n G N and all . . . ,j(n) G N, e(l), . . . , e(n) G {1,*}, corresponding mixed 
moment is uniformly bounded, i.e. I^dl^i ^'(ij)l — a ™ -^ or some non-negative real a n ; 

(4) </? factors over the naturally ordered products in in the sense of Definition^ 

Assume that for all i 7^ j and all e,e' G {1, *}, bf" and b 6 ^ satisfy the commutation relation 



Ve>,e(j,i)bj b[, 



(10) 



Theorem 2 (Extended Non-commutative CLT). Consider a noncommutative probability space 
(A, if) and a sequence of elements {&i}j g N in A satisfying Condition [s| Fix q G M., t > and 
let {fi(i,j)}i<i<j be drawn from a collection of independent, identically distributed, non-vanishing 
random variables, with 

W.iF^ei, E(fx(i,j) 2 ) = i. (11) 

Letting fJL*,*(i,j) = fJ>(i,j) for 1 < i < j, populate the remaining /j, et€ i(i,j), for e,e' G {1,*} and 
i^j € Nj, &y {j) and ©. 
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Then, for almost every sequence j)}i<j , the following holds: for every n G N and all 
e(l),...,e(2n)€{l,*}, 

lim <p(sfK..sf n -V) = 0, (12) 



lim ipiSp . . . Sf n) ) = y q cross(t) t nestCn TT ^(to*)^*)) ^ 

•/ / e,# 2 (2n) 



with Sn G A as given in ([!]), ~f = {(w\, z%), . . . , (w n , z n )}, and where crossi^f) denotes the number 
of crossings in Y and nesti^f") the number of nestings in "V (cf. Definition^. 

The generalized commutation structure of Theorem [2j compared to Theorem [TJ generates a 
second combinatorial statistic in the limiting moments - that of nestings in pair partitions. This 
refinement also extends to the Fock-space level, with the above limits realized as moments of 
creation and annihilation operators on the (q,t)-Fock space [3], briefly overviewed in Section [2j 
Compared to ([5]), the operators of interest now satisfy the commutation relation 

a q ,t(f)a q:t (g)* - q a q>t (g)* a q>t (f) = (/, g) x t N , (14) 



where N is the number operator. Note that (11), together with the fact that t > 0, recovers the 
fundamental constraint that \q\ < t in order for the (q, i)-Fock space to be a bona fide Hilbert 
space. 

The extended non-commutative CLT of Theorem [2] requires that ip{bibi) = ip(b*b*) = (p(b*bi) = 
(cf. Condition [2]), not generally needed in Speicher's setting (Condition [IJ except in the asymptotic 
models for the (/-commutation relati on (Jo"| ). Remark [T] of Section [3] and Remark [4] of Section [4] discuss 
the existence of the limits ( |l2"| ) and (|13[) and the form they take in the absence of this requirement. 
Nevertheless, this additional condition is in fact consistent with the natural choice of matrix models 
- specifically, extending Lemma [TJ the following is a generalized Jordan- Wigner construction. The 
underlying probability space (A n ,^Pn) remains that of the previous section. 

Lemma 2 (Two-parameter Jordan- Wigner Transform) . Fixq G R, t > and let {fi ete /(i, j)}i^j^,e'&{\,*} 
be a sequence of commutation coefficients, i.e. a sequence of non-zero real numbers satisfying Q 
and Q. Consider the 2x2 matrices {(7 X } X ^, 7 given by 





" 1 




"01" 




Vix 


7 = 






For i = 1, ... ,n, let ^{i,j) ■= fJ-*,*(i,j) and consider the element b n ^ G M.%^)® n given by 

K,i = a fi(i,i) ® o>(2,i) <8> ... <8> <V(i-i,j) <8> 7 ® gi g ■ ■ ■ g <7i • (15) 

®(n — i) 



T/ien, /or ewery n £ N, the non-commutative probablity space (A n ,(p n ) and the elements b nj i,b ni 2, 
■ ■ ■ , b n ,n G A n satisfy Condition^ 

Finally, analogously to p], Theorem [2] and Lemma [2] together yield an asymptotic random matrix 
models for the creation, annihilation, and field operators on the (q, t)-Fock space: 

Corollary 1. Consider a sequence of commutation coefficients drawn according to Theorem^ and 
the corresponding matrix construction of Lemma [I| Let 

. Nk 

$N,k '■= —7= ^2 b]S[k,i- (16) 
V i=N(k-l)+l 

Then, for any choice of k, i(l), . . . , i(k) G N, e(l), . . . , e(k) G {1, *}, 

lim VNk (S £ V . . . sfU = Vq 4a qtt (eiY {l) . . . a q , t (e k )^) (17) 

iv— )-oo iv/ iv/ 



and 



JV-£c VM^WM + ^(1)) • • ' ( S N,i(k) + S *N,i(k))) = Pff,t( a «,t(ei) ■ ■ • s qA e k)) (18) 

where <f> q ± is the vacuum expectation state on the (q,t)-Fock space ^ q ^{J^f), operator a q ^{e.i) is the 
twisted annihilation operator on ^ q ,t(J^f) associated with the element of the orthonormal basis 
of Jtf, and s 9i t(ej) := a qi t(fii) + a q ,t{ei)* is the corresponding field operator. 

2. Preliminaries 

The present section overviews the key combinatorial constructs used to encode the mechanics of 
the non-commutative Central Limit Theorem. It also overviews the Hilbert space framework that 
provides a natural setting in which to realize the limits of the random matrix models of Corollary [TJ 

2.1. Partitions. Denote by &(n) the collection of partitions of [n] := {1, . . . ,n}. Set partitions 
will be extensively used to encode equivalence classes of products of random variables, based on the 
repetition patterns of individual elements. Specifically, any two r-vectors will be declared equivalent 
if element repetitions occur at same locations in both vectors; i.e. for (i(l), . . . , i(r)), (j(l), • • • ,j( r )) £ 
[N} r , 

(»(l),...,i(r))~ (j(l),...,i(r)) for all 1 < h < k 2 < r, 

i(k 1 ) = i(k 2 )iSj(k 1 )=j(k 2 ). (19) 

It then immediately follows that the equivalence classes of "~" can be identified with the set &(r) 
of the partitions of [r]. An example is shown in Figure [I] Note that writing "(i(l), . . . , i(r)) ~ Y v 
will indicate that (i(l), . . . , i{r)) is in the equivalence class identified with the partition Y G &(r). 



1 1 m i °° 1 1 m i 

( 8, 8, 1, 8, 1, 6 ) ( 1, 1, 3, 1, 3, 10) 

Figure 1. Two elements of [N] r (for N = 10, r = 6) that belong to the same 
equivalence class, where the latter is represented as the corresponding partition 
r G &{r) given by V = {(1, 2, 4), (3, 5), (6)}. 

Particularly relevant is the collection ^2(2^) of pair partitions of [2n], also referred to as 
pairings or perfect matchings, which are partitions whose each part contains exactly two ele- 
ments. A pair partition will be represented as a list of ordered pairs, that is, <^2(2re) 3 V = 
{(wi, Zi), . . . , (w n , z n )}, where W{ < Zi for i E [n] and w\ < . . . < w n . In the present setting, the 
pair partitions will typically appear with additional refinements given by the following two statistics 
on ^ 2 (2n). 

Definition 1 (Crossings and Nestings). For Y = {(wi, z\), . . . , (w n , z n )} G ^2(2^), pairs (wi, z,j) 
and (vjj,Zj) are said to cross if Wi < Wj < z.- L < Zj. The corresponding crossing is encoded by 
(wi,Wj,Zi, Zj) with Crosstyf) := Wj, z%, Zj) \ (wi, zi), (wj, Zj) 6 "V with wt < Wj < Zi < zj} as 
the set of all crossings in "V and cross^) := \ Cross{'f)\ counting the crossings in Y . 

For Y = {(wi, z\), . . . , (w n , z n )} S ^2(2n) ; pairs (wi, Zi) and (wj, Zj) are said to nest if wi < 
Wj < zj < Z{. The corresponding nesting is encoded by {w{ ) withNest(Y) := {(wi 

(wi, Zi),(wj, Zj) £ Y with Wi < Wj < Zj < Zi} as the set of all nestings in Y and nest(Y) := 
\Nest{Y)\ counting the nestings in Y . 

The two concepts are illustrated in Figures [2] and [3j by visualizing the pair partitions as collections 
of disjoint chords with end-points labeled (increasing from left to right) by elements in [2n]. 
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Figure 2. An example of a crossing [left] and nesting [right] of a pair partition 
V = {(ei,zi), . . . , (e n ,z n )} of [2n]. 




FIGURE 3. Example of three pair partitions on [2n] = {1,...,6}: cross(^i) = 
3, nest(^i) = [left], cross(y 2 ) = 2, nest(r 2 ) = 1 [middle], cross(r 3 ) = 
0, nest(^j) = 3 [right]. 



2.2. Operators on the (g,t)-Fock space. The (q,t)-Fock space ^ q ^{J>f) [3], for \q\ < t, is a two- 
parameter deformation of the classical Bosonic and Fermionic Fock spaces. Consider the tensor 
algebra on the Hilbert space J$? (taken as real and separable) given by ^(Jt?) := ©^(C®^)®", 
with (C ® Jf?)° defined as a complex vector space spanned by a real unit vector Jrff. The 
algebra ^(Jti?) is spanned by the pure tensors {hi (g) . • • <8> h n \ n € N, hi, . . . , h n € J^f} U {0,}. The 
completion of ^(Jif) with respect to the usual inner product, denoted ( , )o and given by (SI, fi)o 
and (/i . . . /„, /ii . . . h m ) = <5 n ,m(/l, ^l)jr • • • {fn, K),%', yields the full (Boltzmann) Fock space. In 
the present scenario, it will be more interesting to complete with respect to the "((/, i)-symmetrized" 
inner product ( , ) q> t given by (fi, Q) q j = 1 and 



I w, 77 ^ 777 

{fl®...®f nt hi®...®hm) qt t={ Jn) t \ I f h \ ( 20 ) 




71 = 771 



where inv(-7r) denotes the number of inversions of the permutation ir £ S n . The completion of 
J^~(^f) with respect to ( , ) q j yields the (q,t)-Fock space ^ q! t(J^), where letting t = 1 recovers the 
g-Fock space ^ q (J4?) of Bozejko and Speicher [JJ. 

The annihilation operators {a q! t(h)}h & ^ on ^q it (Jif ) and their adjoints (with respect to ( , ) ?) t), 
the creation operators {a qt t(h)*}h^jf, are densely defined on by 

a qtt {f)*n = f, a q , t (f)*h 1 ®...®h n = f®h 1 ®...®h n (21) 

n 

o,, t (/)n = 0, a q>t (f)hi ® ...®h n = J2 q k ~H n - k (f, h k )^ hi ® . . . ® h h ® . . . ® h n , (22) 

fc=i 

where the superscript hk indicates that h^ has been deleted from the product. Letting t be 
the linear operator defined by t N Q = $7 and t hi <8> ■ • • <8> h n = t n h\ ® . . . ® h n , the creation 
and annihilation operators are readily shown to satisfy the (q, ^-commutation relation (14). The 
two-parameter family of the (self-adjoint) field operators s qt t(h) := a q j(h) + a q j(h)* , for h £ Jf, 
is referred to as a (q, t)-Gaussian family. The moments of the creation, annihilation, and field 
operators are computed with respect to the vacuum expectation state ip q j ■ £$(^ qt t(J4?)) — > C, 
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ip q j(a) = ( a Q, Q) q ,t- In particular, for every n £ N and all e(l), . . . , e(2n) £ {1, *}, 

p q ,t(a q AhiT {1) ■ ■ ■ a q ,t(h2n-i) e{2n - 1] ) = 0, (23) 

n 

<MMM e(1) • • • a^(^2n) e(2n) ) = Y, g cross(r) t nest(r) J] p(M^J^a,, t (/»*)<*<>), 

(24) 

^(^(/i) 2 "- 1 ) = 0, (25) 
tp q ,t(s q , t (h) 2n ) = \\h\\% J2 g cross(r) i nest(r) , (26) 

where ^2(2^) is again the collection of pair partitions of [2n] and each "V £ ^2(2^) is (uniquely) 
written as a collection of pairs {(wi, zi), . . . , {w n , z n )} with ioi < . . . < w n and Wi < Z{. 

3. The Extended Non-commutative Central Limit Theorem 

The goal of this section is to extend the "deterministic formulation" of the non-commutative 
Central Limit Theorem of Speicher [9]. The deterministic result differs from the previously stated 
Theorem [T] in that the sequence of commutation signs (s(i,j))ij, taking values in {—1,1} and 
associated with the commutation relations 6f6^- = s(j,i)bjbl, is now fixed. In [9], an analogous 
Wick-type formula is nevertheless shown to exist, provided the existence of the following limit: 

j}^ II s(i( Wj ),i(w k )) :=Xy 

i(l),...,i(2n)e[A r ] s.t. (wj ,w k ,Zj ,z fc )eCross(>') 
(i(l),...,i(2n))~r 

for each pair partition ~f £ £P(2n). 

At present, the focus is on a sequence (6j)j 6 N of non-commutative random variables satisfying a 
more general type of commutation relations, where for all i 7^ j and e, e' £ {1, *}, 

b\b £ - = He'Mi^Vjbi for some /v, e (t,i)€R. (27) 

At the outset, the sequence of commutation coefficients {f^e,e'(i, j)}i^j,e,e'€{i,*} m ust satisfy certain 
properties. In particular, interchanging the roles of i and j in the commutation relation implies 
that 

fJ>e,e>{i,j) = (A) 

Similarly, conjugating (via the * operator) both sides of the commutation relation yields 

(For example, 6j6j = fiiA(j,i)bjbi and therefore 6*6* = (6^)* = fii y i(j,i)b*b*, but also 6*6* = 
V*,*(i,j)b*b*.) Therefore, by (A), 

1 1 

= t^v, M*,i(^i) = T^v- ( B ) 

The second key ingredient in a non-commutative CLT is a moment-factoring assumption. As in 
[9], the factoring is assumed to follow the underlying partition structure. Drawing on the notation 
of Section [2J viz. the equivalence relation "~" on the set [N] r of r-tuples in [N] := {1, . . . , N}, the 
two relevant ways in which the moments may be assumed to factor are defined as follows. 

Definition 2. Consider a sequence {bi}i£fq of random variables, elements of some non- commutative 

probability space (A, if). The element 6^*j . . . 61™%, for e(l), . . . , e(n) £ {1, *} and i(l), ■ ■ ■ , i{n) £ N, 

is said to be an interval-ordered product if , . . . , i(n)) ~ Y where "f = {{1, . . . , kx}, {ki + 
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1, . . . , A^}, • • • , + 1 . . . , fci^i}} is an interval partition of [n]. The same element is said to be 

a naturally ordered product if, in addition, i{l) < i(k\ + 1) . . . < i{k\-y\-\ + 1)- 

The state ip is said to factor over naturally (resp. interval) ordered products in {6j}jgN if 

(tfW b e(nh_ (b e(l) e(k^_ 1 + l) e(k m ) 

■■■°i( n )) ■■■ i(k 1 ))---^°i(k lr} _ 1 +l) ••• i(k\-y\) } 

whenever . . . is a naturally (resp. interval) ordered product. 



The following remark ensures that the commutation relations (27) are consistent with the moment 
factoring assumptions. 

Remark 1. In assuming ip factors over naturally ordered products, one must be able to bring a 

moment tp(b^b^b^bj) for i > j into naturally- ordered form. Alternatively, should it be further 
assumed that ip factors over interval- ordered products of the sequence {bi}i^fq, one must allow that 

(p(b^ l b i i bj i bj') = <p(bj bj b^bf) for all i,j and e, e' € {1,*}- When commutation coefficients are 

e' e e ' 

constrained to take values in { — 1, 1}, it is in fact the case that b^b^ commutes with b-b-, and the 
moment-factoring assumptions are consistent with the commutativity structure. However, this need 
not be the case for the general setting. In particular, 

m^bjbf ) = //,.,■ (j, /;//, ,, (j, /)//,<.,; (j, /)//,;,. (j, i) <p(bpfb?b<). 

The reader may verify that any sequence of real-valued commutation coefficients for which the above 
product evaluates to unity regardless of the choice of e, e' must in fact take values in {—1, 1}. 

Instead, rather than imposing additional restrictions on the sign sequence, the alternative ap- 
proach is that of restricting the range of ip when applied to the sequence {h}. In particular, by 
(A)-(B), 

Me 1 ,e5(i>*)Me il e j O'>OA*<,ejC7>0/*<,ej(i>*) = 1 

whenever e 7^ e' . Thus, by imposing that (p(b*b*) = ip{bibi) = for all i £ N, the assumption 
on the factoring of naturally- ordered second moments conveniently becomes equivalent to factoring 
of interval-ordered second moments. Note that factoring an interval-ordered product containing 
higher moments generally still incurs a product of commutation coefficients. However, as will 
become apparent shortly, the contribution of such expressions vanishes in the limits of interest. 

The stage is now set for the main result of this section. 

Theorem 3 (Extended Non-commutative CLT). Consider a non- commutative probability space 
(A, tp) and a sequence of elements of A satisfying Condition [s| with the real-valued com- 

mutation coefficients {fi e >^(i,j)} satisfying the consistency conditions (A)-(B). For n £ N, fix 
e(l), . . . , e(2n) G {1,*} and, letting ^2(2n) denote the collection of pair partitions of[2n], assume 
that for all Y = {{w\, zi), . . . , (w n , z n )} 6 <^2(2ra) the following limit exists: 

( 



V,e(l),...,e(2n) : = J im N " Y] 



i(l),...,i(2n)e[N] s.t. 
(i(l),...,i(2n))~f 



II ^e{ Zj ),e{w k ){i{Zj)^{ w k)) 



(Wj,Wk,Zj,Z k ) 

\ &Cross{V) 



\ 



II ^( 2J ), e ( Zm )( i ( 2; i)^( 2: m))Me( 2i )^(«' m )( i ( Z i)' i ( U, ™)) 

(wj,Wm,Z m ,Zj) 

G Nest(T) ) 



(28) 



where Cross(Y) and Nest(Y) denote, respectively, the sets of crossings and nestings in Y (cf. 
Definition^ and where the equivalence relation ~ is given by (19). 



Then, for every n G N and all e(l), . . . , e(2n) G {1, *}, 

lim ^...^ 2 ^) = 0, (29) 

iv — >oo 

n 

\^{Sf ] ...Sf n) ) = J2 V, e( l),..., e (2n)n^( 6£K) ^ l) )' ( 30 ) 

re# 2 (2n) i=l 

/or S'at G A as given in and with each "V G ^2(2n) written as ~f = {(wi, Z\), . . . , (w n , z n )} for 
wi < . . . < w n and Wi < Z{ (i = 1, . . . , n). 

Proof of Theorem [3]. The notation and the development follow closely those of [U] • 
Fix r G N and e(l), . . . , e(r) G {1,*} and recall that the focus is the ./V — > oo limit of the 
corresponding mixed moment of Sn. Namely, let 

M N := as** . . . Sp) = Yl <P OS • • • 4i) ■ 

i(l),...,i(r)£[N] 

Making use of the previously-defined equivalence relation, can be rewritten as 

M »= Y. jk E K^'S •••<$) = E 



where 



f£#(r) i(l),...,i(r)e[N] s.t. fe#(r) 

(i(l),...,*(r))~f 



i(l),...,i(r)e[A r ] s.t. 
(i(l),...,i(r))~r 

Focusing on M'^, suppose first that ~f contains a singleton, i.e. a single-element part {A:} G V 
for some k G [r]. Via the commutation relation (10), b^l . . . b^} can be brought into a naturally 
ordered form (incurring, in the process, a multiplying factor given by the corresponding product 
of the commutation coefficients). In turn, by the assumption on the factoring of the naturally 

ordered products (cf. Definition j^J) , ip ^(lj ■ ■ -fr^)) factors according to the blocks in V. Since 

f(bk) = <p(bk) = 0, it follows that for all JV£N, Mjj = for all partitions Y containing a singleton 
block. 

Focus next on partitions of [r] containing blocks with two or more elements or, equivalently, 
partitions "V G ^(r) with \V\ < [r/2\, where \Y\ denotes the number of blocks in "f. Recalling 
that, by the assumption on the existence of uniform bounds on the moments, we have that for all 
f G 0»(r), 



< a. 



for some ay G M. Thus, for a partition 'f with t blocks, summing over all • • • , i(r) G [N] with 
(i(l), . . .,i{r)) ~ f yields 



\M%\ < r £ \£\ ay, 
and therefore 

^ ( N ) /I 

E ^£<*- 

Noting that (1) the above sum is taken over a fixed (finite) index r, (2) that the only TV-dependent 
term in the above expression is the ratio (^)/iV r / 2 and (3) that (^)/-V //2 - > as N — > oo for 
£ < \r/2] , it follows that only those partitions Y with \y\ > [Y/2] contribute to the N — > oo limit 
of Afjy. But, since |^| < L r /2j> it follows that the only non-vanishing contributions are obtained 
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N 



for r even and partitions with exactly r/2 blocks — i.e. pair-partitions, 'f G £? > 2(r). Therefore, for 
r odd, 

lim v(sfK..Sp) = 

N— >oo 

and, otherwise, 



lim \M N \ = Y lim AT" I V p ( 6 e f!j . . . bf\ 

■f6# 2 (r) i(l),...,i(r)6[AT] s.t. 

(i(l),...,i(r))~f 



Next, fixing i(l), . . • , z(r) £ [iV] and recalling that ^ is a pair-partition of [r], consider the following 

e(l) e(r) \ I 

algorithm for transforming b}J . . . b},(, via the commutation relation (10), into an interval-ordered 
product. Starting with and recalling that "V is a pair-partition of [r], let 1 < k\ < r denote the 
unique index for which i{ki) = Consider element b^ to be already in place and commute 

frj(fci) with the elements to its left until b^i) is immediately to the right of recording all the 
while the commutation coefficients incurred in each transposition. The next iteration, proceeding 
in the analogous manner, is carried out on the string of length r — 2 given by i (2), ... , i(k\), . . . , i(r), 
where i{k\) indicates that i(k\) has been suppressed from the string. Continuing in this manner, 
the algorithm terminates when the remaining string is the empty string. The resulting moment is 
of the form 

( b e(l) b e(r)\ = d e(l),...,e(r) ( b e(w!) b e( Zl ) e(w r/2 )e(z r/2 )\ 
^ ■ • • °i(r) ) ~ Pi(l),...,i(r) V \° l (w 1 )°i{z 1 ) ■ ■ ■ °i(w T/2 )°i(z r/2 ) ) ' 

where ~f = {(wi, Z\), . . . , (w r / 2 , z r /2)} with w\ < . . . < w r / 2 is the underlying pair-partition and 

e(l) e(r) 

^i(l) i(V) denotes the product of the commutation coefficients incurred in this transformation. 
Note that though i(wj) = i(zj) for all j = 1, . . . , r/2, in general e(wj) ^ and the above 

expression therefore also (artificially) distinguishes between i{wj) and i(zj). 
While it need not be the case that i{w\) < . . . < i(w r /2), and the moment 



(V(«<l)^Ol) Xu>r/2) b <Z r /2) 

r \ i( Wl ) i(zi) • • • i{w r/2 ) i(z r/2 ) 



therefore need not be naturally ordered, ip nevertheless factors over the pairs. Specifically, as 
ip(bjbj) = ip(b*b*) = 0, it can be assumed that e(wj) 7^ t(zj) for j = 1, . . . , r/2. By Remark [TJ it 
then follows that 

*P \ i(l) ■ ■ ■ °i(r) ) - Pi{l),...,i{r) V { i( Wl fi( Zl ) J ■■■ < P {°i(w r/2 fi(z r/2 ) 

Next, P^m'"''t{r) can be expressed combinatorially as follows. Fixing some j G [r/2] and con- 
sidering the corresponding pair (wj,Zj) G "V (where Wj < Zj), note that for every k G [r/2] for 
which Wj < Wk < Zj < Zk, the above algorithm commutes Zj and Wk- Additionally note that this 
commutation is performed exactly once, on the j th iteration, as the process does not revisit pairs 
that were brought into the desired form in one of the previous steps. The corresponding contribu- 
tion to is therefore given by ii^ Zj )^ Wk ){i(zj),i(w k )). Similarly, for every m G [r/2] for 
which Wj < w m < z m < Zj, the above algorithm commutes Zj and z m as well as Zj and w m , and 

e(l) e(r) 

both commutations occur exactly once. The corresponding contribution to PJ^' "i( r ) ^ s therefore 
given by Me( Zj -),e(2 m ) i{z m ))ii e ^ ie ^ m ^[i{zj),i{w m )). Recall now (cf. Definition [T| that the 
4-tuple given by Wj < w k < Zj < z k is what is referred to as a crossing in Y and encoded by 
(wj,Wk,Zj,Zk) G Cross('^ / ), whereas the 4-tuple Wj < w m < z m < Zj is referred to as a nesting in 
Y and is encoded as (wj,w m , z m , Zj) G Nest('^ / ). Finally, realizing that the algorithm performs no 
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previous iterations commute 



Figure 4. The process of bringing a mixed moment into a naturally-ordered form 
involves commuting all the inversions and all the nestings in each of the underlying 
pair partitions. In commuting a crossing (wj,Wk,Zj,Zk), as depicted, the corre- 
sponding moment incurs a factor ^ e (zj),e(w k ){^{ z j)-,K w k))- 




Naturally ordered on 
previous iterations 



o o 



o - o doobob o o 



Unaffected 



o- o Sh o dTo^b o o 



Unaffected 



Figure 5 . The process of bringing a mixed moment into a naturally-ordered form 
involves commuting all the inversions and all the nestings in each of the underlying 
pair partitions. In commuting a nesting (wj,w m , z m , zj), as depicted, the corre- 
sponding moment incurs a factor M e (^), e ( Zm ) K^))Pe(zj),e(w m ){K z o)ii( w ™))- 



other commutations than the two types described, it follows that 

o6(l),...,e(r; 
\l),...,i(r) 



Pint-Mr) = II ^(z 3 ),e( Wk ){i{Zj),i{w k ))^ 



{wj ,w k ,zj ,z k )eCross(f) 

]J Me(*j),e(^)(*(«i)>K^))Me(* j ),e(io m )(«(«j)>*( 1lw m))- 

)eNest(r) 

The encoding of ^xj' ' '^) through nestings and crossings of "V is illustrated in Figures ^ and J^J 
Putting it all together, 

lim ud(^ {1) Q< 2n h - V lim AT§ V lo ( h e[wi) h e{zi) \ to ( h e{Wn) h e(z "A fl 6 ( 1 )>-' e ( r ) 

N™oo^ [ N N >~ ^ N™oc ^ ^% 1 ) 6 ! ( 2l )j'--^ 6 ! W (, ! ( 2 „)j%..., ! (r) 

re^ 2 (2n) i(l),...,i(r)e[N] s.t. 

(i(l),...,i(r))~y 

(31) 

By the assumption on the covariances of the bi's and the existence of the limit in (128]), 



lim ip(Sf } . . . Sf n) ) = V if ( b< w ^b<^) ...ip( lim AT-I V 

N— too — ' V /V / TV— »oo — ^ ^ /v! v / 

re.3*2(2n) i(l)....,i(r)e[N] s.t. 

(i(l),...,i(r))~f 



which yields ( 30 ) and completes the proof. □ 

Remark 2. The assumption of Theorem^ that the covariances of the bi's are independent of i, 
namely, that 99(6^ 6- 2 ) = (p{b e ^b e ?) for all i,j £ N and e%, €2 £ {1, *}, was not used in obtaining (29) 
and (31). Provided the existence of the limit in (31), the additional assumption is solely used for 
the purpose of simplifying (3l\ ) as (30). 

The above Theorem [3] differs from Theorem 1 of [9] in the following ways: 
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(1) The more general commutation relation b\b e - = /v )e (j,i)b € -b\ with fj, e>e /(i,j) G M now re- 
places the commutation relation bjbj = s(i,j)b € -b\ with spins s(i,j) G {—1, 1}- 

(2) For the purpose of factoring naturally ordered second moments as interval-ordered second 
moments, it is presently additionally assumed that ip(b*b*) = ip{bibi) = 0. (cf. Remark [Tj) 

(3) The convergence of the moments now hinges on the existence of a more complicated limit, 
which is not only a function of the commutation coefficients and of the underlying partition, 
as was the case in [3], but also on the pattern of adjoints in the mixed moment of interest 
(i.e. on the string e(l), . . . , e(n)). 

Note that the assumption on the uniform bounds on the moments is not new, but is instead implicit 
in [9j. 

4. Stochastic Interpolation 
Recall that, analogously to Theorem 1 in [9], the "deterministic version" of the non-commutative 



CLT hinges on an existence of the limit (28), which is determined by the sequence of commutation 
coefficients {fj, e e '(hj)}- Rather than providing more explicit conditions for the existence of the 
above limit, this section follows the philosophy of [PJ and instead considers the scenario where 
the coefficients "may have been chosen at random". The outcome will be that, starting with 
a probability law for a single coefficient and extending it to a product measure on the entire 
coefficient sequence, almost any choice of commutation coefficients will yield a finite and easily 
describable limit. For this, it is first necessary to define a suitable product measure on the coefficient 
sequence that is consistent with the dependency structure given by (A)-(B), which is accomplished 
in Remark[3j In turn, Remark [4] considers the effect on the limit achieved by imposing the vanishing 
of certain second moments. Finally, Lemma [3] is the remaining ingredient in the "almost sure 
version" of the non-commutative CLT (viz. the present Theorem [2]). 

Remark 3. Defining a measure on the sequence of commutation coefficients by focusing on the 
triangular sequences j)}i<i<j o,nd attempting to fix the remaining coefficients via (A)-(B) 

still leaves one degree of freedom. Namely, fJ>*,*(i,j) was not until now explicitly related to j). 
The need for a third relation governing the sign sequence comes into play when considering positivity 
requirements. Generally, <p is assumed to be positive, that is, if ip(aa*) > for all a G .e/. Then, 
<p(bib*) > and ip(bibjb*b*) > for all i, j G N. But, by the commutation relations and the factoring 
of naturally ordered moments, 

ip(bibjbjb*) = f^,i{i,j)fi*,*(i,j)(p(bibi)(p(bjbj). 

If the sequence &i, 62, . . . is such that (p(bib*) > for all i, the commutation signs must therefore 
also satisfy the following, third, requirement: 



I fJ-*,l{hJ) I I V*,*{h3) I 
In the random setting, (C) translates to fJ>*i(i,j) = l(i, j)fJ-*,*(i, j) for some random sequence 
ilihj)} supported on (0,oo). 

In assuming {"f(i,j)} to be independent of fi>**{i,j) in line with the general philosophy of this 
section, the reader may soon verify that only the expectation of 7(1, j) will matter from the perspec- 
tive of Lemma^ Furthermore, since the expectations o/yu* ! *(i,j) and /^*,i(i,j) will be taken to not 
depend on the index (i,j), one is free to set t := Then, for i < j, (C) becomes: 

M*,l(*>i) = *M*,*(»>j)) 1 > 0- ( c ') 



Remark 4. Beyond the existence of the limit (28), the goal of the present section is to develop a 
probabilistic framework in which this limit takes on a particularly natural form. For this purpose, 
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the basic setting of Theorem [3| will need to fulfill an additional requirement. Specifically, by the 
assumption of factoring of naturally- ordered moments, ip(bibjb*b*) and (p(bib*b*bj) for i < j are 
both brought into their naturally ordered form by performing a single commutation. In the former 
case, the commutation incurs a factor [J,*,i(i,j) and, in the latter, the factor Yet, in the 

combinatorial formulation, both products are in the equivalence class ( in the sense of "~ ") of the 
pair partition ir = \ l — i — 1 and both are brought into their naturally ordered form by commuting the 
single crossing in tt. Thus, in order for the combinatorial invariant to be preserved, either: 

• the expected values of j) and /U* 5 *(i,j) must be the same, 
or, 

• one of the two mixed moments vanishes, i.e. ip(b{b*) = or (p(b*bi) = for all i G N. 

By Remark^ one may without loss of generality let = t [i*^{i,j). Thus, as the reader 

may soon be able to verify, opting to make equal the means of j) and ^* t *(i,j) by letting 

t = 1 reduces the statistics of the desired limit to those of crossings and the outcome is the same 
as in the case of randomly chosen commutation signs in |5] . The formulation of Lemma [3] instead 
opts for the second alternative, and the introduction of the second parameter t will give rise to the 
appearance of a second combinatorial statistic, that of nestings. 

Note that while (p is assumed to be positive, it is not assumed to be faithful, and there is no 
contradiction in assuming that (p(b*bi) = while <p(bib*) ^ 0. As further discussed in the following 
section, letting (p(b*bi) = and ip{bib*) = 1 will provide an asymptotic model for a family of 
"twisted" annihilation operators, whereas making the opposite choice would yield the corresponding 
analogue for the creation operators. 

Lemma 3. Fix < \q\ < t and let j)}i<i<j be a collection of independent, identically dis- 
tributed non-vanishing random variables, with 

E( f i(i,j)) = qt- 1 £1, E(>(i,i) 2 ) = l. 

Letting fj.*^(i,j) = n(i,j) for 1 < i < j, populate the remaining (J- e ,e'(hj) f or £) e ' 6 by 

Atl,l(«,j) = T^T, VlAhJ 



I ■ -\ ) V") J I I ■ -\1 



Then, for any ~f £ ,^(2^) and e(l), . . . , e(2n) € {1, *}, the limit \ 28) exists a.s. Moreover, if V 
is such as to satisfy (e(w),e(z)) = (1,*) for all blocks (w,z) G Y, the corresponding limit is given 
by 

V,e(l),...,e(2n) = Q 't a.S. 

where cross(Y) = \ Cross{'V)\ and nest(^V) = \Nest{'f)\ denote, respectively, the numbers of cross- 
ings and nestings in y (cf. Definition^. 

Proof. Fix "V = {(u>i, z\), . . . , (w n , z n )} and, recalling that i(w m ) = i(z m ) for all m 6 [n], consider 
the (classical) random variable Xn given by 

X N := N~ n Yj II ^e{z 3 ),e{w k ){i{Zj),i{Wk)) 

j(l),...,j(2n)e[A r ] s.t. \(wj,w k ,Zj,z k )eCmss(f) 
(i(l),...,i(2n))~r 

X II ^(^), e ( Zm )(^(^) ) i ( 2: m))Me(^),e( ? « m )(^^)^(^m)) , (32) 

(w e ,w m ,Zm,zi)&Nest(i / ) j 

where the sequence of random variables {[i e ^>(i, j)} e ,e'^{i,*},i,j(3^i^j is obtained by letting ^*.*(i, j) = 
fj.{i,j) and fixing the remaining coefficients as prescribed by Q-Q. The first goal is to compute 
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E(Xjv). By the independence assumption, since the overall product includes no repeated terms, the 
expectation factors over the products. It therefore suffices to evaluate ^(fJ- e (z),e(w k )(*(' u; i)' K w k))) 
for each crossing (wj,w k ,Zj,z k ) and E(p e (z t ),e(z m )( i ( w e), ^m))/^)^^) W^), for each 

nesting (wg,w m , z m , Z() of a given pair-partition. At the outset, recall that every pair-partition 
Y contributing to is such that (e(w) , e(z)) = (1,*). Then, starting with the crossings and 
assuming that i(wj) = i(zj) < i(w k ) = i(z k ), the corresponding commutation coefficient in (32) 
and its expectation are given as 

fj,* t l(i(zj),i(w k )) = t^^(i(zj),i(w k )) = t/j,(i(zj),i(w k )) i — > t(q/t) = q. (33) 
whose expectation is t(q/t) = q. 

When it is instead the case that i(wj) = i(zj) > i(w k ) = i(z k ), it suffices to notice that by (A)-(B), 
fJ>*,i(hj) = M*,i(i>0- The same conclusion then holds and each crossing therefore contributes a 
factor of q on average. Moving on to nestings, let (wg, w m , z m , zg) be a nesting. If i(wg) = zg < 
i(w m ) = i(zm), the corresponding commutation coefficient in (32) and its expectation are given as 

^*,*{i(zg),i(z m ))ix %1 (i(zg),i(w m )^ je^ - . 

= t(fi(i(wi),i(w m ))) 2 ' { ' 

If on the other hand i(wg) = Z£ > i(w m ) = i(z m ), by (A)-(B) the commutation coefficient and its 
expectation become 

li*,*(^zi),i(z m ))y>*,i(i(zi),i(w m )) = (fJ-*,*(i(zm),i(z e ))y 1 ^ :1 (i(zi),i(w m )) je^ , . 

= (pi(i(z m ),i(zt)))~ l tii(i(z m ),i(zi)) =t 

Thus, each nesting also contributes a factor of t. It follows that E(Xjv) is given by 

E(X N ) = N- n g cross ( r ) t Dest( - T) = q cross Cn t ^st(T) N -n n , ^ 



i(l),...,i(2n)e[AT] s.t. 
(i(l),...,i(2n))~r 



Thus, 



lim E(X N ) = g cross W f«tW, (37) 

It now remains to show that lirnzv->oo -^iV = limjv_ >00 E(X^) a.s., that is, that for every r] > 0, 
^ frhv>l Uj\/>at{I^m — E(Xm)| > = 0- The calculation is analogous to that in [9j. By the 
subadditivity of P and a standard application of Markov inequality, 



P ( |J {\X M - E(X M )\ >r,}\< F (\ X M ~ E(X M )| > r,) < \ ^ E(\X M - E(X 

\M>N J M>N ^ M>N 

In turn, E(\X M - E(X A /)| 2 ) = H X m) ~ ^(X M ) 2 , with E(X M ) 2 = q 2cmss ( r ) t 2nest ^ and 



m)| 2 )- 
(38) 



E(Xlf) = M~ 2 ™ ^ E 



i(l)....,i(2n)e[AT] s.t. 

(i(l),...,i(2n))~TC, 
j(l),...,i(2n)£[JV] s.t. 

0*(l),-,j(2n))~r 



[| fM*,i(i(^fe)>*(^))) (^*,tU( w k),j(wi)) 



(wk,wa)£ Cross ("^) 



(39) 



J| ^t,t(*(«'m),i(Wn))/'* : lWWm),«(w„))) U** ,* 0' (^m) , j K) ) C? ( w m) , j (w n ) ) 

(u) m ,ttJ m )eNest(T) 

where, for convenience of notation, each crossing (w k ,wg, z k , zg) was abbreviated as (w k ,wg), and 
similarly for the nestings. Now suppose that for two choices of indices and the corresponding sets 
(not multisets) {i(l), • ■ ■ ,i(2n)} and {j(l), ■ • • ,j(2n)}, there is at most one index in common, i.e. 
suppose that {i(l), . . . ,i(2n)} D {j(l), . . . , j(2n)} < 1. In that case, (i(k),i(k')) ^ (j (m) , j (m!)) 
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for all k,k' ,m,m' G [2n] with k ^ k',m ^ m! . By the independence assumption, the above 
expectation factors over the product (up to the parenthesized terms) and the contribution of each 
such . . . ,i(2n)},{j(l), . . . ,j(2n)} is simply q 2cmss ("n t 2nest (' rs ) . Thus, the choices of indices 

with {i(l), . . . ,i(2n)}n{j'(l), . . . ,j(2n)} < 1 do not contribute to the variance E(\X M - E(X M )\ 2 )- 



It now remains to consider the Q(M 2n 2 ) remaining terms of the sum (39). 

By the Cauchy-Schwarz inequality, the expectation of the product is bounded, and thus 

E(|X M - E{X M )\ 2 ) < M~ 2n M 2n ~ 2 C = ^, 

where C does not depend on M. Since Y1m>o M~ 2 converges, 

lim V E(\X M -E(X M )\ 2 ) ->0, 

TV— >oo * — » 

M>N 



and therefore by ( 38 ) 



f| \J {\X M -E(X M )\> V } = HmP (J {|X M -E(X M )|>r/} = 0. 

V -/V>1 M>N J \A/>JV / 

This completes the proof. □ 

5. Random Matrix Models 

Considering some prescribed sequence {fJL e>e r(i, j)} e ,e'E{l,*},i,jeN,v£j °f real- valued commmutation 
coefficients satisfying (A)-(B), Lemma [T] exhibits a set of elements of a matrix algebra that satisfy 
the corresponding commutativity structure. The construction is analogous to the one given in 
Lemma [T] and the latter is in fact stated in a form that renders the present generalization natural. 

Proof of Lemma^ To show that b e - = fJ, e ',e{ji i) bj bj, it suffices to consider the definitions in (15) 
and commute 2x2 matrices. Specifically, let i < j and, by elementary manipulations on tensor 
products, note that 

(<fy(i,i) <g> ... <g> cr,j,(i-i,i) <8> 7 ® vi <8) ... <8> cri (8) erf ) x 

® ••• ® ® ® ® ••• ® 7 ® crf {n ~ j) ) 

= ® • • ■ ® (7^(i,i)) ® (^(i+lj)) 8> • • ■ ® (^17) ® (o-icri)®( n ^ (40) 

Now note that o^o^ = a y a x for all x, y G R. Moreover, 70"^. = \fix a x j. Thus, 
7°>(i,i) = V^M(«) 3)°n{ij)l and 017 = (v / t)~ 1 7o-i, 



and, therefore, 



Next, in commuting b* with bj, the only non-trivial commutations are that of 7* with cWij) and 
o~\ = o~\ with 7. Since ^*a x = (Vix) _1 ovy*, it follows that 

£>*&,• = — = -p - — -bjb* = — ,. . x bjb*A = .. .. bM = ,, .. bM = a\JjA)bM. 

The remaining relations now follow by taking adjoints, and the result is that bfbj = /v,e(j\ i) bj'bf. 

It remains to show that, in addition to the commutation relation, the resulting matrix sequences 
also satisfy the assumptions (l)-(4) of Theorem [2] Start by noting that for a\,...,ak G ^#2, 
<p(ai (8) ... (8) cifc) = (ai)n . . . (ajfc)n, where (a)n := (eia, ei)2- It therefore immediately follows 
that for all i € N, = fib*) = 0. By the same token, it is also clear that for all i,j G N, 
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^(bib*) = (fibjb*) = 1 and y»(6|6f ) = <^(6|6|') = for e,e' £ {1,*} with (e,e') / (1,*), and, 

furthermore, l^dliLi — ^ ^ or an n ano - an choices of exponents and indices. The factoring 

over naturally ordered products also follows immediately, completing the proof. □ 

Finally, combining Theorem [2] with Lemma [2j and comparing the resulting moments with those 



given in Section 2.2 immediately yields the desired asymptotic models for the creation, annihilation, 
and field operators on the (q, i)-Fock space. For instance, the mixed moments of Sn converge to 
those of the annihilation operator a(ei), where e\ is an element of the orthonormal basis of Jtf. 



More generally, the expressions of Section 2.2 in fact consider systems of operators, e.g. they specify 
the joint mixed moments of annihilation operators a(ei), . . . , a(e n ) associated with basis elements 
ei, . . . , e n . In order to asymptotically realize the joint moments of a(ei), . . . , a(e n ) rather than the 
moments of a(e\) alone, it suffices to consider a sequence Sjv,i> • • ■ > SW,n of partial sums built from 
non-intersecting subsets of {&j}j g N- For instance, the fact that e« and ej are orthogonal for i ^ j 
and that the moment <p q j(a(ei)a(ej)) vanishes follows (in this asymptotic setting) from the fact 
that (p{bibj) = for i ^ j. The general formulation is found in Corollary 1. 
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